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1  INTRODUCTION 


In  this  note  we  address  some  questions  related  to  acoustic  coupling  be¬ 
tween  transmitting  elements  in  an  acoustic  array.  We  shall  be  primarily 
concerned  with  elements  of  characteristic  size  a  that  are  small  compared 
with  the  acoustic  wavelength  A: 


a«A.  (1-1) 

We  shall  also  suppose  that  the  array  is  operated  in  a  linear  acoustic  regime. 
This  will  be  interpreted  as  implying  that  the  acoustic  pressure  p(a)  at  the 
transducer  is  less  than  pcav,  that  at  which  cavitation  occurs: 

p(a)<Pcav-  (1-2) 

The  array  is  assumed  to  have  NT  radiators  with  phase  centers  at  locations 

X{(i  =  1  ,...Nt)  in  the  neighborhood  of  the  origin  of  a  rectangular  coordi¬ 
nate  system.  Taking  advantage  of  condition  Equation  (1-1)  we  compromise 
slightly  on  generality  by  assuming  spherical  radiators  of  radius  a  that  are 
driven  in  a  spherically  symmetric  mode.  Then  the  acoustic  pressure  at  a 
distance  R,  from  the  ith  element  is 

Pl(r)=p,(a)(a/Rt)e^R'-a\  (1-3) 


where 

R,  =  r  -  x, . 

Since  we  have  assumed  a  linear  system,  the  total  acoustic  pressure  at  r  is 

P(r)  =  Y^P.ir).  (1-1) 

i=i 

In  the  far  field  of  the  array  we  can  write  this  as 

P(r)  =  p„(a/R)c'kRD(r),  (1-5) 

where,  to  within  an  unimportant  phase  factor, 

D(r  )=£()-, (1-6) 

1=1 

! 


and,  say, 

P.  =  £  |p,(n)hV,.  (1-7) 

J 


I  he  maximum  radiated  beam  power  will  occur  in  a  given  direction  r  when 
the  phases  of  the  p,(a)  are  adjusted  so  that 


D(  r)  =  Nr. 

In  this  case  the  maximum  radiated  power/unit  solid  angle  is 

Pbeam  =  r2\p\2/{2pc) 

-  °-^N2 

~  2  pclR- 

[Here  p  =;  1000  kg/m3  and  c  1500  m/s  is  the  speed  of  sound.] 

The  total  radiated  power  is 

p  _  /  l/)(r)[2  2  j0 

/fad  J  2 pc  r  dn' 

where  the  integral  is  over  all  solid  angles  SI  .  Wc  obtain 

p  .  _ 


(1-8) 


rad 


W /v>i  r  , 


(1-9) 


(1-10) 


E.,  ,  ,  .  sin[A|x,  —  x,  1]  . 

pAaMa)  l|x,  -x,|  )• 

Of  special  interest  is  a  linear  or  planar  array  for  which  the  radiated  beam 
is  "broadside,”  or  r  •  x,  =  0.  Then,  if  maximum  beam  power  is  desired  the 
driving  phases  will  be  equal,  so 


-Itra2  f;Vr 

^  rad  =  El/>.(«)|2 

t  i=i 


(Ml) 


E,  .  ...  sin  [A:  lx,  -  x.j] 


When  high  beam  power  is  the  objective,  the  second  term  in  E 
(  1-10)  is  not  of  great  significance  (see  Section  2)  and  we  have 


ion 


/ 


rad 


l7Trt2  „  , 


<  l-n2.Yr 


/Pcav\ 

V  ~PC  )  ' 


(M2) 


For  purposes  of  illustration,  we  shall  take 


PcaV/(2pc)  =  13/v  W/m\  (1-13) 

Prad  <  4tui2  NTx[l3I\W/m2}.  (1  -  14) 

Returning  to  Equation  (1-8),  we  conclude  that 

/’beam  <  «* N?x\nKW /m1].  (1-15) 


The  bounds  in  Equations  (1-12)  and  (1-15)  on  the  total  radiated  power 
[see  the  clarification  in  Section  2]  and  on  the  beam  power /steradi an  appear  to 
be  fundamental  for  simple  arrays  operating  in  a  linear  regime.  The  quadratic 
dependence  on  element  radius  indicates  the  price  which  is  paid  in  total  radi¬ 
ated  energy  in  exchange  for  the  price  saved  with  small  radiators. 

Our  argument  docs  not  bound  radiators  working  in  a  nonlinear  regime  or 
radiators  placed  near  reflectors,  baffles,  etc. 

In  Section  2  we  examine  the  role  of  the  second  term  in  Equations  (1-10) 
and  (1-11).  In  Section  3  we  investigate  the  performance  of  an  array  when 
close  acoustic  coupling  between  radiating  elements  obtains. 


:l 


2  THE  TOTAL  RADIATED  POWER 


In  this  Section  we  consider  the  effect  of  element  displacement  in  changing 
the  total  radiated  power  in  an  array.  For  a  linear  or  planar  array  with 
broadside  steering  we  can  replace  Equation  (1-11)  by 


P' 


'  +  f£ 
T 


sin[A;jxi  —  x;  | 
fc|x,  -X_,| 


(2-1) 


dropping  multiplicative  factors. 


First  we  consider  a  linear  array  with  A/2  element  spacing.  Then  P'  =  1. 
Evidently,  we  can  maximize  Equation  (2-1)  by  setting  all  the  x,  equal.  Then 
P'  =  Nr,  but  we  have  no  beam  forming!  Various  more  restrictive  algorithms 
have  been  tried,  none  of  which  increase  beam  power  [consistent  with  Equation 
(1-8)].  We  illustrate  this  for  a  10  element  array  and  a  simple  algorithm  that 
holds  fixed  the  first  and  last  elements  at  positions  0  and  4.5  (we  shall  scale 
all  lengths  in  units  of  A/2).  The  remaining  8  elements  are  sequentially  moved 
between  neighbors  to  maximize  P\  but  constrained  so  no  two  elements  can 
be  closer  than  1/8.  Starting  with  an  element  spacing  of  1/2  and  P'  =  1,  the 
algorithm  converged  in  a  few  iterations  to  the  spacing  shown  in  Figure  1  and 

P'  =  2.57.  (2  -  2) 

The  resulting  beam  pattern  is  shown  in  Figure  2.  The  extra  radiated  energy 
is  fed  into  sidelobes,  as  might  have  been  anticipated  from  Equation  (1-8). 


As  a  second  example  we  consider  a  square  array  of  5  rows,  each  having 
5  elements.  The  spacing  of  rows  and  elements  is  1/2  (i.e.  A/2).  For  this  25 
element  array, 

P'  =  1.0 1 .  (2-3) 

With  sequential  moving  of  elements  between  neighbors  and  rows  between 
rows  (holding  boundaries  fixed  and  permitting  no  elements  or  rows  to  be 
closer  than  1/8)  we  are  led  to  the  element  positions  shown  in  Figure  3  and 

P'  =  1.45.  (2-1) 


Again,  this  extra  radiated  energy  has  gone  into  sidelobes  rather  than 
beam  power. 
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Element  Positions 


Mean  Spacing  =  1/2  Lambda 


Element  Number 


Figure  1  Element  spacing  for  a  10  element  linear  array  to  maximize  the  radiated  power  (2.1) 
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Radiated  Intensity  (db) 
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Figure  2  Beam  pattern  for  the  array  snown  in  Figure  (1). 


3  ACOUSTIC  COUPLING  WITHIN  AN  AR¬ 
RAY 


In  this  Section  we  look  in  some  detail  at  the  beam  forming  mechanism 
for  our  simple  model  of  Nr  identical  spherical  radiators  in  a  linear  or  planar 
array.  The  radius  of  a  radiator  is  a  -f  £,  where 

(  =  w,  the  surface  velocity  (3  —  1) 

and 

ICI  <  a. 

We  also  assume  simple  harmonic  radiation: 

w  =  V  e1^.  (3  -  2) 

'Then  the  acoustic  pressure  at  the  radiator  surface  is 

p{a)  =  Zw  (3  -  3) 

where 

Z  —  p  cka /(ka  -f  ?)  ( 3  —  4 ) 

is  the  acoustic  impedance.  In  obtaining  Equation  (3-3)  we  have  used  the 
relation 

dp /dr  =  iujpw 


and  Equation  ( 1-3). 


f  or  t  he  array  of  A  r  elements  the  average  pressure  at  the  surface  of  element 


i  is 


where 


,Vr 


J= 1 


"b  =  VJ ( 


(3  -  5) 

(3  —  t>) 


I  he  average  pressure  in  Equation  (3-5)  is  defmed  as 

,)=  [PP 

i  i-..1 


m-ps 
1  • 


9 


(3-7) 


where  the  integration  extends  over  the  spherical  element  surface.  Using 
Kquations  ( 1  -3), (3-3),  and  (3-7)  we  obtain 


y  y  « 

=  /.-T< 
" IJ 


iK»„-n)si  »(*•«) 
ka 


(3-8) 


We  note  at  this  point,  that,  the  relation  Equation  (3-5)  neglects  scattering 
of  radiation  l>v  the  radiators.  This  scattering  is  dependent  on  the  specific 
array  geometry  and  is  expected  to  be  small  for  small  radiators.  For  the 
sphere  i  and  radiation  from  j  we  have  [see  M.  J linger  and  D.  Feit,  Sound. 
Structures  and  Their  Interactions.  MIT  Press,  1986] 

^cattered  %  {ka)^a/ RJ.  (3  -  9) 

'  radiated 

At  any  rate,  we  shall  neglect  sound  scattering  within  the  array. 

1  lit'  equations  for  driving  the  array  surface  are  assumed  to  have  the  form 

A/ ic,  +  a  ic,  =  Ft  -/>,(«).  (3-10) 


Here  A/  is  a  mass  and  o  is  a  damping  constant  characterizing  the  radiators, 
and 

F,  =  f,  c-‘- 

On  substituting  Equation  (3-5)  into  Equation  (3-10)  and  removing  the  time 
dependent,  factors,  we  obtain 


[-'■->  A/  +  a  +  X\\\  =  /,  -  £  (;A  -  1  1 ) 

.)(#<) 


Now,  the  hn(Z)  represents  amass  of  water  and  the  Rt(z)  radiation  damping, 
so  we  may  write 


-In, (7.) 
II' (Z) 


=  A/„  = 


1  Tva'p 


(ka  )2  + 

(ka)M,l. 


1  hus.  we  ran  re-write  Ecpiation  (3-1  1)  as 

[  '-A/,  +»']W  =/,  - 


CM3) 


(3  -  13) 


10 


Here 


Me  —  M  -j-  Mw 
ae  =  a+u(ka)Mw 

=  wir? i +  (*«)"»] 

and  vve  have  introduced  the  of  the  radiator: 

Mu 


Q  = 


4  7m 


The  pressure  at  great  distances  from  the  array  is 

Nr 


p{r)  =  Z  \  z)  T,exp[ik(Ri  -  a)]K- 

1=1 


(3-14) 


(3-15) 


(3-16) 


It  is  convenient  to  write 


/,  =  |/i|e'a'  (3-17) 

Vi  =  \Vi\eu-. 

To  steer  a  beam  in  the  horizontal  direction  (from  broadside)  0,  we  set 


c,  =  A\r,sin<?3. 

The  power  input  to  the  array  is  (for  T  <^u~l) 

1  rT 


Pin  =  fJo  {£/>,}<« 

=  Hl^ll/|cos(a,  -  e,)/2 


(3-18) 


(3-19) 


and  the  total  power  radiated  is 


1=1 


(3-20) 


The  radiated  power  in  the  beam/unit  solid  angle  is  obtained  using  Equations 
(1-5)  and  (3-16) 

r2\p(r)\2 


I  beam 


2  pc 


(3-21) 


!  1 


For  an  array  for  which  acoustic  coupling  can  be  neglected  we  use  Equation 
(3-18)  and  set  |l/,j  =  V’u.  a '  1  i.  Then  from  Equation  (3-13)  we  find 

/,  =  +  o,] i;  (3-22) 

These  driving  forces,  witii  the  phases  defined  by  Equation  (3-18),  will  be  used 
as  "standard”  in  Equation  (3-13)  when  we  are  interested  in  the  coupling. 

Eor  numerical  illustration  we  shall  consider  two  arrays.  The  first  is  an 
array  of  3  rows  each  having  7  elements  in  a  horizontal  line.  The  second  array 
has  5  rows  each  having  5  elements  in  a  horizontal  line.  Half  wavelength 
spacing  is  assumed  in  both  arrays.  We  shall  further  take 

Q  =  0.2,  Mw  =  M,  and  uj  =  628.  (3  -  23) 

With  Equation  (3-22)  as  the  specified  driving  force,  we  can  solve  the  coupled 
equations,  Equation  (3-13)  for  tin*  element  surface  velocities  V,  .  The  far 
field  pressure  is  finally  obtained  from  Equation  (3-16).  We  emphasize  that 
for  all  calculations  shown  here  the  driving  forces  are  obtained  from  Equation 
(3-22)  with  the  phases  set  by  Equation  (3-18). 

We  first  consider  radiat  ion  from  the  7x3  element  array  when  the  elements 
are  uncoupled  [  i.e.,  ZtJ  =  0  for  i  ^  j  in  Equation  (3-23)].  For  broadside 
steering  (i.e.  0,  =  0)  the  radiated  intensity  as  a  funct  ion  of  observation  angle 
0  (in  the  horizontal  plane  and  measured  from  0S  is  shown  in  Figure  T  The 
intensity  as  shown  is 

HO)  =  ioiog(|/;|2/A;2),  (3  -  2i) 

where  l)  is  given  by  (1-6).  In  Figures  -5  through  8  we  show  the  radiated 
intensity  obtained  from  coupled  elements  using  Equation  (3-13).  These  are 
shown  for  different  radiator  scaled  sizes: 

dia  =  la/ A.  (3  —  2o) 

In  Figures  !!  and  10  we  show  the  radiated  intensity  for  0S  =  30°. 

In  the  next  series  of  figures,  we  show  the  corresponding  radiated  power 
from  the  5  x  5  array. 

In  looking  at  these  figures,  we  note  that  increasing  the  coupling  strength 
(i.e.,  the  element  radius  a)  reduces  tin1  effective  acoustic  impedance  and 


Received 
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Acoustic  Array  Intensity  7x3  Elements:  Uncoupled 


Figure  4  The  radiated  power  (3.24)  for  the  7x3  array  as  a  function  of  horizontal  angle  when  the 
steering  angle  0S=  0. 


Radiated 


Acoustic  Array  Intensity  7x3  Elements  Coupled 


Figure  5  The  7x3  array  radiated  power  for  coupled  elements.  9S=0,  and  an  element  diameter  of  0  38. 
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Acoustic  Array  Intensity  7x3  Elements:  Coupled 
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Figure  6  The  radiated  power,  as  in  Figure  (5),  but  for  an  element  diameter  0.64. 


Radiated 


Acoustic  Array  Intensity  7x3  Elements:  Coupled 


30  -15  0  15  30 

Observation  Angle  (degrees) 


Figure  7  The  radiated  power,  as  in  Figure  (5),  but  for  an  element  diameter  0.89. 


Hi 


Radiated  Intensity  (db) 
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30  -15  0  1  5 

Observation  Angle  (degrees) 


Figure  8  The  radiated  power,  as  in  Figure  (5),  but  for  an  element  diameter  0.96. 


Received 


Acoustic  Array  Intensity  7x3  Elements:  Uncoupled 


Figure  9  The  radiated  power,  as  in  Figure  (4),  but  for  0S  =30°. 
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Acoustic  Array  Intensity  7x3  Elements:  Coupled 


Figure  10  The  radiated  power,  as  in  Figure  (5),  but  for  0S=  30\ 
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Acoustic  Array  Intensity  5x5  Elements:  Uncoupled 
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Figure  1 1  The  radiated  power  (3.24)  for  the  5x5  array  as  a  function  of  horizontal  angle  when  the 
steering  angle  0S=0  and  the  elements  are  uncoupled. 
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Acoustic  Array  Intensity  *>*.5  Elements  Coupled 


Figure  13  The  radiated  power,  as  in  Figure  (12),  but  for  an  element  diameter  0.64. 


Radiated 


Acoustic  Array  Intensity  5x5  Elements:  Coupled 


Figure  14  The  radiated  power,  as  in  Figure  (12),  but  for  an  element  diameter  0  89. 
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Figure  16  The  radiated  power,  as  in  Figure  (11),  but  for  9S  =  30°. 
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Acoustic  Array  Intensity  5x5  Elements:  Coupled 


Observation  Angle  (degrees) 


Figure  17  The  radiated  power,  as  in  Figure  (12),  but  for  0S  =30°. 
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Radiated  Intensity  (db) 


Acoustic  Array  Intensity  5x5  Elements:  Coupled 


Figure  18  The  radiated  power,  as  in  Figure  (15),  but  for  0S  =  30° 


Radiated 


Acoustic  Array  Intensity  5x5  Elements  Coupled 


Figure  19  The  largest  element  excitation  (3.26),  normalized  to  that  for  the  uncoupled  system, 
is  shown  as  a  function  of  element  diameter. 


Beam  Power(coup).  Beam  Power(uncoup) 


Diameter 

Figure  20  The  beam  power  is  shown  as  a  function  of  element  diameter. 


2!) 


increases  the  radiated  power  (holding  fixed  the  driving  forces).  This  is  illus¬ 
trated  in  f  igure  1!)  where  the  largest  element  excitation 


|\ |2max 

Hol2 


(3-26) 


is  shown  as  a  function  of  the  element,  diameter.  The  corresponding  beam 
power  is  shown  in  Figure  20.  Perhaps  the  most  obvious  effect  of  beam  cou¬ 
pling  is  to  reduce  the  minimum  radiated  power  between  beam  lobes-a  not 
surprising  result,  as  this  results  from  a  delicate  balance  of  phases. 


For  a  linear  array,  the  effects  of  element  coupling  can  in  principle  be  com¬ 
pensated  for  by  adjusting  the  driving  forces  /,.  The  practical  consequences 
of  coupling  for  beamforming  do  not  seem  to  be  very  great,  however. 
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